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We perform a global parton fit to DIS and related data, including next-to-leading 
logarithmic (NLL) BFKL resummations in both the massless and massive sectors. The 
resummed fit improves over a standard next-to-leading order (NLO) DGLAP fit, with 
a positive definite ghion at the input scale as opposed to the negative gluon seen at 
NLO. Furthermore, the predicted longitudinal structure function is free of perturbative 
instability at small x, and the reduced cross-section shows a turnover at high y (absent 
in the NLO fit) consistent with the HERA data. 

1 Small x Resummation 

1.1 Motivation 

Current and forthcoming particle collider experiments involve very high energies, such that 
the momentum fractions x of initial state partons are extremely small. The splitting func- 
tions that govern the evolution of parton densities fi{x,Q 2 ) with momentum scale Q 2 , to- 
gether with the coefficients that relate these partons to proton structure functions, are unsta- 
ble at low Bjorken x values due to terms behaving like x _1 a§ log m (l/a;) where n > m + 1. 
Although the standard DGLAP theory (where the splitting and coefficient functions are 
considered at a fixed order in as) works well in parton fits, there is some evidence that a 
resummation of small x logarithms is necessary. Previous work has shown that a LL analysis 
fails to describe data well. One resums small x logarithms in the gluon density by solving 
the BFKL equation [2], an integral equation for the unintcgrated gluon 4-point function. 
One then relates this gluon to structure functions using the kr factorisation formalism [3l [4] 
to obtain the resummed splitting and coefficient functions. 

1.2 Solution of the BFKL equation 

Introducing the double Mcllin transformed unintcgrated gluon density: 




(1) 



the NLL BFKL equation in (JV, 7) space is a double differential equation in 7: 



d-y 2 



d 2 fi(j,Q 2 ) 
d~i 2 



1 rf(xo( 7 )/(7,JV)) 
P N d 7 
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with Po = 3/(7r/3o). The derivatives in 7 arise from the use of the LO running coupling 
as(k 2 ) = l/(/?ologfc 2 /A 2 ) in momentum space, and Xn{l) is the Mellin transform of the 
n th -order BFKL kernel. One may solve this to give: 

-"(-1$) W)* (2) 

for some A(j) and Xi("f). One would ideally like to factorise the perturbative from the non- 
perturbative physics to make contact with the collincar factorisation framework. This can 
be achieved (up to power-suppressed corrections) by shifting the lower limit of the integral 
in equation (|2|) from 7 — > 0. Then one finds for the integrated gluon: 

G(N,t) =G E {N,t)Q I {Ql,N), (3) 

where the perturbative piece is: 

1 /.1/2+JOO rf3 

Gh(N, t) = — — exp [jt - Xi( 7j N)/(poN)] d 7 , (4) 

where X\ can be derived from xo(7) and Xi(7)> an d is a known function of 7. Structure 
functions have a similar form: 

Fh(N,t) = — l7 ' j/ cxp[ 7 f-X 1 (7,iV)/(/3 jV)]d 7l (5) 

2 ™ Jl/2-too 7 

where /i( 7 , TV) is a NLL order impact factor coupling the virtual photon with the BFKL 
gluon. If all impact factors are known, one can derive all necessary splitting and coefficient 
functions in double Mellin space (within a particular factorisation scheme) by taking ratios 
of the above quantities. The non-perturbative dependence then cancels, and one obtains 
results in momentum and x space by performing the inverse Mellin integrals either numeri- 
cally or analytically. The exact NLL impact factors are not in fact known, but the LL results 
supplemented with the correct kinematic behaviour of the gluon have been calculated [51 [3] . 
We have shown that one expects them to approximate well the missing NLL information in 
the true impact factors [7]. 

Consistent implementation of small x resummations in the massive sector requires the def- 
inition of a variable flavour number scheme that allows the massive impact factors to be 
disentangled in terms of heavy coefficient functions and matrix elements. We have devised 
such a scheme, the DIS(x) scheme jSj. With resummations in both the massive and mass- 
less sectors, one has everything necessary to carry out a global fit to DIS and related data. 
First, the resummed splitting and coefficient functions are combined with the NLO DGLAP 
results using the prescription: 



plot. _ pl\ LL p 



NLL 1 u NLO 



pNLL(O) _|_ pNLL(l) 



where the subtractions remove the double counted terms, namely the LO and NLO (in 
as) parts of the resummed results. Then the resulting improved splitting and coefficient 
functions interpolate between the resummed results at low x, and the conventional DGLAP 
results at high x. 
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2 Results 



The resummcd splitting functions P+ (~ P gg + A/9P gg at small a;) and P qg are shown in 
figure [U One sees that the LL BFKL results are much more divergent than the standard 
NLO results, which are known to describe data well. The addition of the running coupling 
to the LL BFKL equation suppresses this divergence, but it is still unacceptable. Inclusion 
of the NLL BFKL kernel, however, leads to a significant dip of the splitting functions below 
the NLO results. This dip is also observed in other resummation approaches [HI [TU] and 
has an important consequence in the global fit in that it resolves the tension between the 
Tevatron jet data (which favour a larger high x gluon) and the HI and ZEUS data (which 
prefer a higher low x gluon). By momentum conservation, one cannot increase the gluon 
at both low and high i in a standard NLO DGLAP fit. This is possible in the resummed 
approach, due to the dip in the splitting functions. 




Indeed, the gluon distribution at 
the parton input scale of Q\ = 
lGeV 2 is positive definite over the 
entire x range. This is in con- 
trast to a NLO fit, where the gluon 
distribution is negative at small x 
for low Q 2 values. Whilst a nega- 
tive gluon is not disallowed, it can 
lead to negative structure func- 
tions which are unphysical. The 
resummcd gluon, however, leads to 
a prediction for the longitudinal 
structure function that is positive 
and growing at small x and Q 2 , 
in contrast to fixed order results 
which show a significant pcrturba- 
tive instability. 

A consequence of a more sensi- 
ble description for Fl is that a 
turnover is observed in the reduced 
cross-section a — F 2 — y 2 /[l + (1 — 
y) 2 ] Fl at high y. As seen in fig- 
ure [51 this is required by the HERA data. Furthermore, this feature is missing in NLO fits 
(but present at NNLO). Thus the resummations lead to qualitatively different behaviour, 
consistent with known consequences of higher orders in the fixed order expansion. Overall, 
we find very compelling evidence of the need for BFKL effects in describing proton structure 

HD. 
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Figure 1: Splitting functions in the DIS scheme for 
n/ = 4, t = log(Q 2 /A 2 ) = 6: NLL+NLO (solid); LL 
with running coupling + LO (dashed); LL + LO (dot- 
dashed); NLO (dotted). 
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Figure 2: Reduced cross-section data, compared with both resummcd and fixed order pre- 
dictions. 
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